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Observations with the Interstellar Boundary Explorer (IBEX) have shown enhanced energetic neutral 
atom (ENA) emission from a narrow, circular ribbon likely centered on the direction of the local interstellar 
medium (LISAA) magnetic field. Here, we show that recent determinations of the local interstellar velocity, 
based on interstellar atom measurements with IBEX, are consistent with the interstellar modulation 
of high-energy (tera-electron volts, TeV) cosmic rays and diffusive propagation from supernova sources 
revealed in global anisotropy maps of ground-based high-energy cosmic-ray observatories (AAilagro, 
Asy, and IceCube). Establishing a consistent local interstellar magnetic field direction using IBEX 
ENAs at hundreds to thousands of eV and galactic cosmic rays at tens of TeV has wide-ranging 
implications for the structure of our heliosphere and its interactions with the LISAA, which is 
particularly important at the time when the Voyager spacecraft are leaving our heliosphere. 


T he heliosphere — the region surrounding 
our solar system that is carved out by the 
supersonic solar wind — moves through the 
local interstellar cloud (LIC), which is the part 
of the galactic environment surrounding the 
Sun. The heliosphere deflects much of the harm- 
ful cosmic-ray radiation from the local galactic 
environment, thereby regulating the radiation envi- 
ronment throughout the solar system, which likely 
has important implications for Earth’s atmosphere 
and life (7). The Interstellar Boundary Explorer 
(IBEX) recently provided updated values for the 
velocity vector of the heliosphere through the LIC 
(2) and direction for the local interstellar medium 
(LISM) magnetic field (Table 1 and Fig. 1) from 
the center of the IBEX ribbon of energetic neutral 
atom (ENA) emission (3-9). These results show 
that the interstellar flow, after removing the mo- 
tion of the Sun through the local standard of rest 
(LSR), is nearly perpendicular (87.6 ± 3.0°) to 
the LISM magnetic field. The LSR describes the 
velocity frame in which the mean motion of the 
oldest stars in the Milky Way in the neighborhood 
of the Sun is zero and is the reference frame in 
which cosmic rays assume near uniformity (in ve- 
locity direction). The local interstellar magnetic 
field direction from the IBEX ribbon center is 
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within ~33 ± 20° of the magnetic field direction 
derived by interstellar polarization data from stars 
within 40 pc (10, 11). Here, we show that the 
anisotropy maps of high-energy (TeV) cosmic 
rays likely provide independent confirmation of 
the interstellar magnetic field orientation inferred 
from the IBEX ribbon center. 

The flux of high-energy (TeV) galactic cos- 
mic rays (GCRs) varies as a function of look di- 
rection in the sky. The large-scale structure in the 
TeV GCR sky consists of two broad asymmetries 
with flux variations of -0.2%: a deficit of GCR 
flux at high galactic latitudes (the “loss cone”) 
and an excess of flux in the heliotail (72) direc- 
tion [the “tail-in” excess (73-27)]. Small-scale 
TeV anisotropies (<~10°) in cosmic-ray arrival di- 
rections possibly arise from cosmic-ray propaga- 
tion in a turbulent magnetic field (22). Because 
TeV GCRs have gyroradii of <700 astronomical 
units (AUs) (<0.005 parsec or pc; 1 pc = 3 x 

Fig. 1. LIC velocities and LISM 
magnetic field directions. The 

LIC velocities inferred from IBEX (2) 
in the heliocentric (HC; green vec- 
tor) rest frame and the local stan- 
dard of rest (LSR; black vector) in 
the galaxy. These vectors are shown 
in both the galactic plane (lower 
plane) and the plane containing 
the galactic poles (upper plane) rel- 
ative to the magnetic field of the 
local interstellar medium (LISAA; blue 
lines). The LIC velocity and magnetic 
field are a part of the Loop I super- 
bubble roughly centered on the 
Lower Centaurus Crux subgroup of 
the Scorpius-Centaurus Association. 

The cosmic-ray anisotropy (gray vec- 
tor) indicates average streaming of 
cosmic rays. 


10 18 cm = 200,000 AU) in the LISM, the ob- 
served GCR asymmetries must originate in the 
immediate interstellar environment of the Sun. This 
study extends and builds on increasing knowledge 
of the LISM magnetic field (3, 4, 10, 11, 23, 24) 
and interstellar flow (2, 25-28) to understand the 
origin of the large-scale high-energy (TeV) cosmic- 
ray anisotropy, which has remained uncertain. 

The observed cosmic-ray anisotropy has a num- 
ber of features that may reflect ordering both by 
the entry of cosmic rays into the heliosphere and 
by the interstellar magnetic field. It has been found 
(73) that the tail-in excess weakens above 10 TeV 
because of the relative sizes of the GCR gyroradius 
and heliotail. Data from underground detectors 
in both hemispheres has led to precise positions 
of the asymmetries (14). Additionally, the under- 
ground muon observatory Super-Kamiokande I 
obtained similar asymmetries for the 10-TeV GCR 
flux, with a deficit toward the constellation of 
Virgo representing the northern galactic deficit, 
and an excess in the direction of the heliotail, 
toward the constellation of Taums (75). The direc- 
tion of the Taums excess is within 29° of the local 
interstellar magnetic field defined by the IBEX 
ribbon (Table 1) and within 23° from the local 
interstellar downwind direction (2). Milagro ob- 
servations in the Northern Hemisphere (19), in 
combination with IceCube observations in the 
Southern Hemisphere (20), have confirmed that 
the galactic deficit is centered at high galactic 
latitudes and extends into both hemispheres (27). 
A global GCR anisotropy model constructed with 
unidirectional and bidirectional components and 
applied to results of the Tibet Asy experiment at- 
tributed the bidirectional flow to GCRs drifting 
parallel to the local interstellar magnetic field. 

The cosmic-ray anisotropy, according to stan- 
dard diffusion theory [e.g., (29, 30)], is proportional 
to the average streaming of cosmic rays. Here, 
we take a small [0.3% (37)] ratio of perpendic- 
ular to parallel diffusion, as well as an interstellar 
flow direction to be nearly perpendicular to the 
magnetic field on the basis of IBEX observations. 
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In this parameter regime, cosmic rays are largely 
guided by the interstellar magnetic field because 
both the anisotropic and perpendicular compo- 
nents of the anisotropy (fig. SI) are small com- 
pared to the field-aligned component. 

We consider two scenarios for the formation 
of the cosmic-ray anisotropy. In the first scenario, 
the magnitude of the anisotropy is determined in 
part by spatial gradients in the average cosmic- 
ray density in response to the interstellar flow 
(eq. S6 in supplementary text). For example, an 
outward plasma flow from the center of the Lower 
Centaums Cmx (LCC) cluster (Fig. 1) pushes 
some cosmic rays out of the local interstellar en- 
vironment in which the heliosphere resides. This 
creates a spatial gradient in the cosmic-ray density 
in the direction of the source and flow, which the 
average streaming (gray vector, Fig. 1) of cosmic 
rays tends to oppose. 

The streaming direction of cosmic rays is de- 
termined on relatively small scales (on the scale 
of the cosmic-ray gyroradii, less than 1% of a pc), 
whereas the streaming magnitude is determined 
over much larger scales of tens of pc (32, 33) 


controlled by the scattering (random changes in 
direction) of cosmic rays caused by interactions 
with the turbulent interstellar magnetic field. 

It is important to differentiate (34, 33) be- 
tween the interstellar magnetic field on local scales 
(thousands of AU or -0.005 to 0.01 pc) and on 
larger (parsec) scales owing to the presence of 
turbulence (36-41). Turbulence disrupts steady 
flows by random motion. In the interstellar me- 
dium, turbulent motion causes tangling and com- 
plexity in the structure of the interstellar magnetic 
field. The coherence scale length in the interstel- 
lar medium — the approximate distance along which 
the magnetic field appears relatively ordered — is 
typically 1 to 10 pc (37). Hence, the average 
magnetic field direction at these scales could be 
very different than in the vicinity of the helio- 
sphere. At the large field-to-flow angle (87.6 ± 
3.0°) found by IBEX (excepting the small region 
within 0.3° of exact perpendicularity, where per- 
pendicular diffusion dominates for the 0.3% ratio 
of perpendicular to parallel diffusion taken here), 
the cosmic-ray density gradient becomes large so 
that the projection of the cosmic-ray streaming 


Table 1. LISM velocity and magnetic field direction. 


Magnitude Galactic Galactic Equatorial right Equatorial 
(km/s) longitude (°) latitude (°) ascension (°) declination (°) 

LISM flow velocity* 

(HC frame) 

LISM flow velocity 
(LSR frame*) 

Interstellar magnetic 
field* 

*Heliocentric (HC) rest frame (2). fHeliocentric velocities are converted to the local standard of rest (LSR) using the solar apex 
motion in (54). JField directions and uncertainty inferred from the IBEX highest-energy steps [1.79 and 2.73 keV (55)] in 
which the ribbon maintains coherence and has the largest line-of-sight. 


23.2 ± 0.3 
18.0 ± 0.9 


185.25 ± 0.24 
47.9 ± 2.9 
210.5 ± 2.6 


-12.0 ± 0.5 
23.8 + 2.0 
-57.1 ± 1.0 


78.5 ± 0.6 
267.0 + 3.0 

48.5 ± 1.5 


18.0 ± 0.5 
23.2 ± 3.1 
- 21.2 + 1.6 
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Fig. 2. TeV cosmic-ray anisotropies compared with predictions. Comparison between observed 
(left) and modeled (right) cosmic-ray relative intensities across the sky (J2000 coordinates). Black 
curves show the magnetic equator with a magnetic field direction derived from the center of the IBEX 
ribbon. On the Left, the region below 25°S latitude is the anisotropy map from IceCube with a median 
energy of 20 TeV ( 18 ) and above 20°S latitude is the anisotropy map from ASy with 5-TeV median energy 
(15). Similarly the modeled map (right) at 20 TeV is shown below 25°S latitude and at 5 TeV above 20°S 
latitude. Both portions of the maps are smoothed over 3° to 5°. Labels indicate upwind and downwind 
directions ( 2 ), the current locations of Voyager 1 (VI), and Voyager 2 (V2) directions, and the "upheld" 
and "downfield" directions. Downfield is along the LISM magnetic field determined by IBEX in the 
direction closest to the interstellar velocity, and upheld is in the opposite direction. Plots are in equatorial 
coordinates with 0 hours at the right and increasing longitudes toward the left. 


into the flow plane opposes the interstellar veloc- 
ity. The cosmic-ray gradient in this case depends 
largely on the ratio of perpendicular to parallel 
diffusion. Reduced levels of perpendicular diffu- 
sion cause increased density gradients and thus 
stronger interstellar modulation of cosmic rays. 

The -0.2% high-energy (TeV) cosmic-ray an- 
isotropy suggests a magnitude for the ratio of 
perpendicular to parallel diffusion of only 0.3% 
(fig. SI). This ratio is compatible with recent re- 
sults (31), but an order-of-magnitude smaller than 
the derived -4% ratio of perpendicular to parallel 
diffusion (41) based on cosmic-ray lifetimes in the 
galaxy (39). Simulations of cosmic-ray diffusion 
(42) similarly predicted a ratio of perpendicular 
to parallel diffusion of -2 to 4% at MeV to GeV 
energies, which are lower than the TeV energies 
considered here. The relatively small ratio of per- 
pendicular to parallel diffusion of -0.3% may be a 
result of lower turbulence levels in the LIC com- 
pared to other portions of the galaxy. 

In our second scenario for the cosmic-ray an- 
isotropy, we consider the contributions of super- 
nova remnants (SNRs) to the local cosmic-ray 
gradient (43-48). A statistical model of cosmic 
rays originating from SNR sources based on the 
temporal and spatial distributions of supernova 
sources (43) suggests that cosmic rays are injected 
locally by different sources and time dependently 
diffuse throughout the galaxy. In the case of the 
convective anisotropy (supplementary text, section 
2), the cosmic-ray gradient is created by steady 
cosmic-ray diffusion against the interstellar flow. 
Therefore, the cosmic-ray gradient driven by SNRs 
depends on the time and spatial distribution of 
SNR sources. Because parallel diffusion in the 
LISM is far more rapid than perpendicular diffu- 
sion, cosmic rays have a strong tendency to stream 
along the interstellar magnetic field. Provided 
that the GCR density gradient along the LISM 
magnetic field has the same sign as that of the 
convective gradient and that the cosmic-ray den- 
sity gradient produces a maximum anisotropy 
magnitude similar to that in observations, the an- 
isotropy associated with diffusion from SNR 
sources has a global morphology similar to that 
of the convective anisotropy (supplementary text, 
section 5). That is, the dominance of parallel dif- 
fusion in the LISM can result in a similar global 
morphology for both the convective anisotropy 
and that driven by SNR sources. 

We developed a model of the LISM magnetic 
field that is deflected around the heliosphere (sup- 
plementary text, section 3) and analyzed its influ- 
ence on high-energy (TeV) cosmic-ray anisotropies. 
We constructed a sky map (Fig. 2) of cosmic-ray 
flux as viewed from the Sun, using Monte-Carlo 
calculations of 10 4 individual cosmic-ray trajec- 
tories in this perturbed magnetic field structure. 

General ordering about the magnetic equator 
deduced from the IBEX ribbon is apparent in the 
high-energy (TeV) cosmic-ray observations (Fig. 2, 
left). Some features absent in the simulation results 
can be attributed to the model’s lack of interstellar 
turbulence (22) that should cause small-scale 
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anisotropic structures, local GCR acceleration 
that may result in a more defined excess near 
the heliotail (21), and asymmetric influences on 
the heliosphere by the interstellar magnetic field 
(49-53). The modeled and observed anisotropies 
show some differences (e.g., the region from right 
ascension = -30 to 90° shows weaker modeled 
than observed anisotropies). 

That a similar ordering of the global cosmic- 
ray anisotropy maps is generated simply by 
taking into account IBEX observations of low- 
energy interstellar neutral atoms to deduce the 
interstellar velocity, and the IBEX ribbon to de- 
duce the local interstellar magnetic field direction 
reveals consistency with independent observations 
across 10 orders of magnitude in particle energy 
(keV energies in the IBEX ribbon compared to 
10-TeV cosmic ray energies). Thus, local inter- 
stellar conditions play a key role in ordering very- 
high-energy (TeV) cosmic rays in the immediate 
interstellar environment of the Sun. 

References and Notes 

1. N. A. Schwadron, H. E. Spence, R. Came, Eos 92, 297 
(2011). 

2. D. ]. McComas et ai, Science 336, 1291-1293 (2012). 

3. D. ]. McComas et ai, Science 326, 959-962 (2009). 

4. N. A. Schwadron et ol., Science 326, 966-968 (2009). 

5. H. 0. Funsten et ai, Science 326, 964-966 (2009). 

6. S. A. Fuselier et ai, Science 326, 962-964 (2009). 

7. ]. Heerikhuisen et ol., Astrophys. ]. 708, L126-L130 
( 2010 ). 

8. ]. Heerikhuisen, N. V. Pogorelov, Astrophys. J. 738, 29 
( 2011 ). 

9. N. A. Schwadron, D. ]. McComas, Astrophys. J. 764, 92 
(2013). 

10. P. C. Frisch et ol., Astrophys. J. 760, 106 (2012). 

11. P. C. Frisch, N. A. Schwadron, http://arxiv.org/abs/ 
1310.2922 (2013). 


12. D. ]. McComas, M. A. Dayeh, H. 0. Funsten, G. Livadiotis, 
N. A. Schwadron, Astrophys. ]. 771 , 77 (2013). 

13. K. Nagashima, K. Fujimoto, R. M. Jacklyn, 7. Geophys. Res. 
103 , 17429 (1998). 

14. D. L. Hall et ol., J. Geophys. Res. 104 , 6737 (1999). 

15. G. Guillian et ai, Phys. Rev. D Part. Fields Gravit. Cosmol. 
75 , 062003 (2007). 

16. M. Amenomori et ai, in Proceedings of the 9th Annual 
International Astrophysics Conference, AIP Conference 
Proceedings, J. Le Roux, G. P. Zank, A. ]. Coates, 

V. Florinski, Eds. (American Institute of Physics, College 
Park, MD, 2010), vol. 1302, p. 285. 

17. M. Amenomori; The Tibet ASy Collaboration, Astrophys. Space 
Sci. Trans. 6, 49-52 (2010). 

18. A. A. Abdo et ai, Phys. Rev. Lett. 101 , 221101 (2008). 

19. A. A. Abdo et ai, Astrophys. ]. 698, 2121-2130 (2009). 

20. R. Abbasi et ai, Astrophys. J. 740 , 16 (2011). 

21. P. Desiati, A. Lazarian, Astrophys. ]. 762 , 44 (2013). 

22. G. Giacinti, G. Sigl, Phys. Rev. Lett. 109, 071101 (2012). 

23. R. Lallement et ai, Science 307 , 1447-1449 (2005). 

24. M. Opher, E. C. Stone, T. I. Gombosi, Science 316, 
875-878 (2007). 

25. M. Witte, Astron. Astrophys. 426, 835-844 (2004). 

26. E. Mobius et ai, Science 326, 969-971 (2009). 

27. E. Mobius et oi, Astrophys. ]. Suppl. Ser. 198, 11 (2012). 

28. M. Bzowski et ai, Astrophys. J. Suppl. Ser. 198, 12 
( 2012 ). 

29. L. ]. Gleeson, W. I. Axford, Astrophys. ]. 154, 1011 
(1968). 

30. M. A. Forman, L. ]. Gleeson, Astrophys. Space Sci. 32, 
77-94 (1975). 

31. A. Shalchi, I. Biisching, A. Lazarian, R. Schlickeiser, 
Astrophys. J. 725 , 2117-2127 (2010). 

32. I. Biisching, M. S. Potgieter, Adv. Space Res. 42, 
504-509 (2008). 

33. V. S. Ptuskin, I. V. Moskalenko, F. C. Jones, A. W. Strong, 
V. N. Zirakashvili, Astrophys. ]. 642 , 902-916 (2006). 

34. ]. R. Jokipii, Science 307, 1424-1425 (2005). 

35. ]. R. Jokipii, Science 316, 839-840 (2007). 

36. L. C. Lee, J. R. Jokipii, Astrophys. ]. 206, 735 (1976). 

37. A. H. Minter, S. R. Spangler, Astrophys. ]. 458 , 194 
(1996). 

38. S. Chandrasekhar, G. Munch, Astrophys. ]. 115 , 103 
(1952). 


39. J. R. Jokipii, E. N. Parker, Astrophys. J. 155, 799 (1969). 

40. J. W. Armstrong, B. ]. Rickett, S. R. Spangler, Astrophys. ]. 
443 , 209 (1995). 

41. J. R. Jokipii, in Proceedings of the 2nd Guillermo Haro 
Conference, J. Franco, A. Carraminana, Eds. (Cambridge 
Univ. Press, Cambridge, 1999), p. 70. 

42. J. Giacalone, J. R. Jokipii, Astrophys. ]. 520, 204-214 
(1999). 

43. P. Blasi, E. Amato, 7. Cosmol. Astropart. Phys. 2012, 011 
( 2012 ). 

44. A. D. Erlykin, A. W. Wolfendale, Astropart. Phys. 25, 
183-194 (2006). 

45. I. Biisching, M. Pohl, R. Schlickeiser, Astron. Astrophys. 
377, 1056 (2001). 

46. M. Pohl, D. Eichler, Astrophys. ]. 766, 4 (2013). 

47. I. Busching et ai, Astrophys. ]. 619, 314-326 (2005). 

48. C. Evoli, D. Gaggero, D. Grasso, L. Maccione, Phys. Rev. 
Lett. 108, 211102 (2012). 

49. G. P. Zank, Space Sci. Rev. 89, 413-688 (1999). 

50. R. Ratkiewicz et ai, Astron. Astrophys. 335 , 363 (1998). 

51. N. V. Pogorelov, G. P. Zank, Astrophys. J. 636, 
L161-L164 (2006). 

52. M. Opher, E. C. Stone, P. C. Liewer, Astrophys. J. 640, 
L71-L74 (2006). 

53. N. V. Pogorelov, J. Heerikhuisen, G. P. Zank, 

S. N. Borovikov, Space Sci. Rev. 143, 31-42 (2009). 

54. R. Schonrich, J. Binney, W. Dehnen, Mon. Not. R. 

Astron. Soc. 403 , 1829-1833 (2010). 

55. H. O. Funsten et ai, Astrophys. ]. 776, 30 (2013). 

Acknowledgments: We thank all those who made IBEX 
possible. IBEX is primarily funded by NASA's Explorers Program 
(Contract no. NNG05EC85C). IBEX data are available at 
http://ibex.swri.edu/researchers/publicdata.shtml. IceCube cosmic 
ray data are available from http://icecube.wisc.edu/science/data. 

Supplementary Materials 

www.sciencemag.org/content/343/6174/988/suppl/DCl 
Supplementary Text 
References (56-67) 

21 August 2013; accepted 29 January 2014 
Published online 13 February 2014; 

10. 1126/science. 1245026 


Nanoporous BiV0 4 Photoanodes with 
Dual-Layer Oxygen Evolution Catalysts 
for Solar Water Splitting 

Tae Woo Kim and Kyoung-Shin Choi* 

Bismuth vanadate (BiV0 4 ) has a band structure that is well-suited for potential use as a 
photoanode in solar water splitting, but it suffers from poor electron-hole separation. Here, we 
demonstrate that a nanoporous morphology (specific surface area of 31.8 square meters per gram) 
effectively suppresses bulk carrier recombination without additional doping, manifesting an 
electron-hole separation yield of 0.90 at 1.23 volts (V) versus the reversible hydrogen electrode 
(RHE). We enhanced the propensity for surface-reaching holes to instigate water-splitting chemistry 
by serially applying two different oxygen evolution catalyst (OEC) layers, FeOOH and NiOOH, which 
reduces interface recombination at the BiV0 4 /0EC junction while creating a more favorable 
Helmholtz layer potential drop at the OEC/electrolyte junction. The resulting BiV0 4 /Fe00H/Ni00H 
photoanode achieves a photocurrent density of 2.73 milliamps per square centimenter at a 
potential as low as 0.6 V versus RHE. 


N -type bismuth vanadate (BiV0 4 ) has re- 
cently emerged as a promising photoanode 
for use in water-splitting photoelectro- 
chemical cells because it absorbs a substantial 
portion of the visible spectrum (bandgap energy, 


-2.4 eV) and has a favorable conduction band 
(CB) edge position very near the thermodynamic 
H 2 evolution potential (1,2). However, the solar- 
to-hydrogen (STH) conversion efficiency achieved 
with BiV0 4 to date has been far below what is 


expected because the material suffers from poor 
electron-hole separation yield (c|) se p) (2-6). Previ- 
ous efforts to improve the 4> sep of BiV0 4 mainly 
focused on doping studies, which were intended 
to improve its poor electron transport properties 
( 2 , 6 - 12 ). 

Here, we demonstrate that a high-surface-area, 
nanoporous BiV0 4 electrode composed of par- 
ticles smaller than its hole diffusion length can 
effectively increase (|) sep without additional dop- 
ing. Furthermore, we investigated the effect of an 
oxygen evolution catalyst (OEC) layer on the in- 
terfacial recombination at the BiV0 4 /0EC junc- 
tion, water oxidation kinetics, and the Helmholtz 
layer potential drop at the OEC/electrolyte junction 
using two different OECs, FeOOH and NiOOH. 
Our understanding of the Bi V 0 4 /OEC/electrolyte 
junction resulted in the development of a new 
strategy to serially apply dual layers of OECs that 
can optimize both the BiV0 4 /0EC and the OEC/ 
electrolyte junctions simultaneously, enabling effi- 
cient utilization of surface-reaching holes for solar 
water oxidation. 
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